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I.  INTRODUCTION 


The  flow  field  in  the  base  region  of  a  jet-propelled  tactical  missile  or 
a  rocket-assisted  artillery  projectile  is  complex.  The  propulsive  jet  is 
underexpanded  and  can  strongly  affect  the  base  pressure  distribution  and  the 
afterbody  flow  field.  For  large  jet  to  free  stream  pressure  ratios,  the 
interaction  of  the  exhaust  jet  with  the  external  flow  can  be  large  enough  to 
induce  extensive  afterbody  flow  separation.  This  can,  of  course,  seriously 
affect  the  control  effectiveness  of  any  control  surface  located  in  that 
region.  For  over  a  decade,  experimental  efforts  and  the  component  approach  of 
Korst1  has  provided  valuable  insight  into  this  flow  problem.  The  flow  field 
under  consideration  is  complicated  due  to  strong  viscous/inviscid  interaction 
and  regions  of  flow  separation.  The  component  approach  has  its  serious 
limitations  for  three-dimensional  configurations  and  similar  procedures  are 
extremely  difficult  to  apply  at  transonic  speeds.2 

Recent  advances  in  numerical  algorithms  and  the  advent  of  supercomputers 
have  made  numerical  modeling  of  the  Navier-Stokes  equations  a  reality.  As 
computational  speed  has  sharply  increased,  the  computing  costs  have  sharply 
dropped.  This  trend  in  technology  makes  the  direct  solution  of  Navier-Stokes 
equations  even  more  attractive.  The  Navier-Stokes  computational  technique 
models  the  strong  interactions  involved  between  the  flow  regimes  in  a  fully 
coupled  manner  and  does  not  contain  the  empiricism  found  in  the  component 
approach.  Some  empiricism  does  enter  into  the  Navier-Stokes  solutions  through 
turbulence  modeling  which  is  an  area  of  further  computational  and  experimental 
research. 

Recently,  Sahu,  et  al,3  computed  the  base  region  flow  field  for  a  pro¬ 
jectile  at  transonic  speeds  using  the  thin  layer  Navier-Stokes  equations.  A 
unique  flow  field  segmentation  procedure  was  used  for  complete  numerical  simu¬ 
lation  of  a  projectile  including  the  base  region.  This  technique  was  also 
used  to  calculate  the  effect  of  a  centered  jet  on  the  base  region  flow  field 
at  a  high  transonic  speed. 4  Other  Navier-Stokes  solutions  have  been  obtained 
for  flow  over  afterbodies  with  exhaust  jet.5-8  The  work  of  Deiwert5  and  Fox6 
are  based  on  the  thin-layer  approximation  of  the  Navier-Stokes  equations 
whereas  Wagner7  and  Thomas8  solved  the  full  Navier-Stokes  equations.  The 
two-layer,  algebraic,  Baldwin-Lomax9  turbulence  model  was  used  in  References 
3-7.  Thomas,  et  al8  used  the  two-equation  k-e  and  k-w  turbulence  models  and 
some  differences  in  the  computed  results  were  observed. 

This  report  describes  the  computational  investigation  of  the  effect  of 
exhaust  jet  on  the  base  pressure  and  the  base  region  flow  field  at  supersonic 
speed  using  the  same  numerical  procedure  of  References  3  and  4.  A  unique  flow 
field  segmentation  procedure,  equivalent  to  using  multiple  adjoining  grids,  is 
used  which  preserves  the  sharp  base  corner.  Additionally,  the  grids  in  the 
base  region  were  adapted  to  the  free  shear  layer  as  the  solutions  developed. 
Solutions  have  been  obtained  for  supersonic  flow  over  an  axisymmetric  boattail 
afterbody  where  the  free  stream  Mach  number  is  2.0  and  the  jet  exit  Mach 
number  is  2.5.  Ratios  of  jet  static  pressure  to  free  stream  static  pressure 
are  considered  in  the  range  of  1  to  9  for  a  conical  nozzle  exit  half  angle  of 
20°.  Comparison  of  the  numerical  results  have  been  made  with  available 
experimental  data. 


II. 


COMPUTATIONAL  TECHNIQUE 


The  Azimuthal  Invariant  (or  Generalized  Axi symmetri c)  thin-layer  Navier- 
Stokes  equations  for  general  spatial  coordinates  £,  n,  5  can  be  written  as10 

3xq  +  3?E  +  3?G  +  H  =  Re_13?S  (1) 

where  £  =  £(x,y,z,t)  is  the  longitudinal  coordinate 


n  =  n(y,z,t)  is  the  circumferential  coordinate 


and 


?  =  c(x,y,z,t)  is  the  near  normal  coordinate 
t  =  t  is  the  time 


p 

PU 

pu 

*  1 
>  e  -  J 

puU+£xp 

pv 

• 

pvU+£yp 

pw 

PWU+£zp 

_e  . 

(e+p)U-£tp 

pW 

PUW+CXP 

pvW+CyP 

pwW+czp 

(e+p)W-ctP 


9 


H  =  J' 


0 

0 


pV[Rc(U-£t)  +  R?(W-?t)] 
-PVR(V-nt)  -  p/R) 

0 


0 


s  = 


‘,<sx+5yt4>u5  *  (»/3)(t<vyt,y!i!« 

“('  xt?y+5x)vc  + 

u(c2+C2+C2)w  +  (p/3)(c  u  +e  v  +e  w  k 
vsx  sy  sz'  ?  '  '  /v^x  5  7  5  sz  c'  z 

{ (c^+Cy+C^JC (m/2) (u2+v2+w2) ^ 

+  KPr"1  (y-1)"1 (a2)^]  +  (m/3) 

(?xu+«yv+«zw)(?Xuc+V?+?zwc)} 


The  velocities 


u  =  h  +  V  +  ^yv  +  Czw 

V  =  nt  +  nxu  +  n^v  +  nzw 

w  *  «t  +  ?xu  +  Cy v  +  Czw 

represent  the  contravariant  velocity  components. 

The  Cartesian  velocity  component  (u,  v,  w)  are  nondimensi onal i zed  with 
respect  to  aw  (free  stream  speed  of  sound).  The  density  (p)  is  referenced  to 

pOT  and  total  energy  (e)  to  p^a2  .  The  local  pressure  is  determined  using  the 

equation  of  state, 


p  =  (y  -  1 ) [e  -  0.5p(u2  +  v2  +  w2)] 


(3) 


where  y  is  the  ratio  of  specific  heats. 

In  Equation  (1),  axisymmetric  flow  assumptions  have  been  made  which 

result  in  the  source  term,  H.  The  details  of  how  this  is  obtained  can  be 
found  in  Reference  10  and  are  not  discussed  here.  Equation  (1)  contains  only 
two  spatial  derivatives.  However,  it  retains  all  three  momentum  equations  and 
allows  a  degree  of  generality  over  the  standard  axisymmetric  equations.  In 
particular,  the  ci  rcumferential  velocity  is  not  assumed  to  be  zero  thus 
allowing  computations  for  spinning  projectiles  to  be  accomplished. 


The  numerical  algorithm  used  is  the  Beam-Warming  fully  implicit, 
approximately  factored  finite  difference  scheme.  The  algorithm  can  be  first 
or  second  order  accurate  in  time  and  second  or  fourth  order  accurate  in 
space.  Since  the  interest  is  only  in  the  steady-state  solution.  Equation  (1) 
is  solved  in  a  time  asymptotic  fashion  and  first  order  accurate  time 
differencing  is  used.  The  spatial  accuracy  is  fourth  order.  Details  of  the 
algorithm  are  included  in  References  11-13. 

For  the  computation  of  turbulent  flows  a  turbulence  model  must  be 
supplied.  In  the  present  calculations  a  two-layer  algebraic  eddy  viscosity 
model  by  Baldwin  and  Lomax9  is  used.  In  their  two  layer  model  the  inner 
region  follows  the  Prandtl-Van  Driest  formulation.  Their  outer  formulation 
can  be  used  in  wakes  as  well  as  in  attached  and  separated  boundary  layers.  In 
both  the  inner  and  outer  formulations  the  distribution  of  vorticity  is  used  to 
determine  length  scales  thereby  avoiding  the  necessity  of  finding  the  outer 
edge  of  the  boundary  layer  (or  wake).  The  magnitude  of  the  local  vorticity 
for  the  axi symmetric  formulation  is  given  by 


-  /  (iih2  +  (J1  _  _£J1\ 
'3z  3y' 


3  V  3W\  2 


3w 

3x 


3Z 


(4) 


In  determining  the  outer  length  scale  a  function9 


F(y)  =  yM  [1  -  exp(-y+/A+)] 


(5) 


is  used  where  y+  and  A+  are  the  conventional  boundary  layer  terms.  For  the 
base  flow  (or  wake  flow)  the  exponential  term  of  Equation  (5)  is  set  equal  to 
zero.  In  other  words,  the  Van  Driest  damping  term  is  not  applicable  and  is 
thus  neglected.  The  outer  formulation  also  requires  the  computation  of  the 
Klebanoff  intermittency  function  and  a  velocity  scale  given  by 

udif  =  <u2  +  v2  +  w2Uax  -  <u2  +  v2  +  w2Uin 


Both  of  the  terms  on  the  right-hand  side  of  Equation  (6)  are  evaluated  via  the 
velocity  profiles.  For  wall-bounded  flows,  the  minimum  term  in  is 

usually  zero.  For  wakes,  the  Klebanoff  intermittency  factor  is  determined  by 
measuring  the  distance  from  the  centerline  of  symmetry.  The  algebraic  eddy 
viscosity  model  may  not  be  strictly  valid  for  all  of  the  wake  flow 
situations.  More  realistic  or  complex  turbulence  models  much  be  considered  a 
subject  left  for  future  study. 


III.  METHOD  OF  SOLUTION 
1.  BASE  REGION  FLOW  WITH  JET-OFF 


The  procedure  used  to  compute  the  base  flow  without  jet  for  a  projectile 
configuration  has  been  described  in  Reference  3;  however,  limited  details  will 


be  repeated  here  for  clarity.  The  code  can  compute  the  full  flow  field 
(including  the  base  region)  of  a  projectile.  However,  for  supersonic  flow  the 
forebody  solution  can  be  obtained  efficiently  using  a  space-marching  Para¬ 
bolized  Navier-Stokes  (PNS)  code.  This  technique  is  used  to  provide  upstream 
boundary  condition  (line  EF,  Figure  1)  for  the  computation  of  the  afterbody 
flow  field  containing  the  base  region.  The  afterbody  solution  is  obtained 
using  the  unsteady  or  time-marching  Navier-Stokes  equations.  Figure  1  shows  a 
schematic  illustration  of  the  flow  field  segmentation  used  in  this  study  for 
computational  purposes.  It  shows  the  transformation  of  the  physical  domain 
into  the  computational  domain  and  the  details  of  the  flow  field  segmentation 
procedure  in  both  domains.  This  flow  field  segmentation  procedure  is  equiva¬ 
lent  to  using  multiple  adjoining  grids.  An  important  advantage  of  this  proce¬ 
dure  lies  in  the  preservation  of  the  sharp  corner  at  the  base  and  allows  easy 
blending  of  the  computational  meshes  between  the  regions  ABCD  and  AEFG.  Mo 
approximation  of  the  actual  sharp  corner  at  the  base  is  made  and  is  inherent 
in  the  current  procedure. 

In  Figure  1,  the  cross  hatched  region  represents  the  model.  The  line  BC 
is  the  base  and  the  region  ABCD  is  the  base  region  or  the  wake.  The  line  AB 
is  a  computational  cut  through  the  physical  wake  region  which  acts  as  a 
repetitive  boundary  in  the  computational  domain.  Implicit  integration  is 
carried  out  in  both  5  and  c  directions.  Note  the  presence  of  the  lines  BC 
(base)  and  EF  in  the  computational  domain.  They  both  act  as  boundaries  in  the 
computational  domain  and  special  care  must  be  taken  in  inverting  the  block 
tridi agonal  matrix  in  the  K  direction.  The  details  of  this  procedure  can  be 
found  in  Reference  3  and  are  not  included  here. 

The  no  slip  boundary  condition  for  viscous  flow  is  enforced  by  setting 


U  =  V  =  W  =  0  (7) 


on  the  body  surface  including  the  base.  Along  the  computational  cut  (AB),  the 
flow  variables  above  and  below  the  cut  are  simply  averaged  to  determine  the 
boundary  conditions  on  the  cut.  On  the  centerline  of  the  wake  region,  a 
symmetry  condition  is  imposed  and  free  stream  conditions  are  used  on  the  outer 
boundary. 


2.  BASE  REGION  FLOW  WITH  JET-ON 

The  method  of  solution  for  the  case  with  a  centered  propulsive  jet 

remains  essentially  the  same  as  described  in  Section  III  A.  The  boundary 

conditions  on  the  body  surface,  at  the  cut  and  the  downstream  boundary  also 
remain  the  same  as  previously  described.  Along  the  base  boundary,  the  same 

conditions  described  earlier  are  used.  For  the  nozzle  exit,  boundary 

conditions  are  used  based  on  the  nozzle  exit  Mach  number-,  stagnation  tempera¬ 
ture  and  pressure.  The  velocity  components  are  linearly  interpolated  from  the 
center  line  of  symmetry  to  the  nozzle  height  at  the  exit,  i.e.,  conical  flow 
at  the  jet  exit  has  been  assumed. 


IV.  MODEL  AND  COMPUTATIONAL  GRID 


The  model  geometry  used  in  this  study  comprised  of  a  2  caliber  14°  half¬ 
angle  conical  nose,  6  caliber  cylindrical  mid-section  and  a  1  caliber  8° 
boattail.  The  nozzle  exit  diameter  is  0.6  calibers.  Detailed  experimental 
measurements  for  this  shape  and  the  same  flow  conditions  described  earlier 
have  been  made  by  Agrell,  et  a l.14 

When  computations  over  the  entire  model  are  made,  only  a  limited  number 
of  grid  points  can  be  used  in  the  base  region.  One  way  to  eliminate  this 
restriction  is  to  use  known  data  given  by  experiment  or  otherwise  at  a  station 
upstream  of  the  base  and  then  compute  the  flow  field  in  the  isolated  base 
region  only.  This,  of  course,  allows  a  large  number  of  grid  points  to  be  used 
in  the  base  region  and  can  be  used  to  determine  grid  dependancy  on  the  com¬ 
puted  solution  in  the  base  region.  This  is  ideally  suited  for  the  numerical 
computations  of  base  region  flow  field  at  supersonic  velocities.  Solutions 
can  be  obtained  for  the  forebody  with  the  space-marching  PNS  code.  In  the 
present  study,  the  PNS  code  was  used  to  generate  a  solution  at  a  station  1.5 
calibers  upstream  of  the  base;  and  this  solution  was  then  used  as  an  upstream 
boundary  condition  for  the  computation  of  the  base  region  flow  field  by  the 
unsteady  base  flow  code. 

Figures  2-5  show  the  expanded  view  of  the  computational  grids  in  the  base 
region.  These  grids  are  shown  for  both  upper  and  lower  halves  for  clarity; 
however,  computations  are  made  only  for  the  upper  half  plane  for  axi symmetric 
flow.  Figure  2  shows  the  grid  for  the  jet-off  case  whereas  the  grids  shown  in 
Figures  3-5  are  for  the  jet-on  cases  of  Pj/Pc  =  1,  3  and  9,  respecti vely. 

Additionally,  these  grids  were  adapted  to  the  free  shear  layer  as  the  solu¬ 
tions  developed.  Logic  has  been  implemented  to  adjust  the  grid  cut  AB  (Figure 
1)  to  the  viscous  shear  layer  which  begins  to  neck-down  shortly  behind  the 
base.  The  height  of  the  cut  is  weighted  between  a  moment  of  shear  and  the 

standard  nonadaptive  grids.  Specifically,  the  cut  height  z,  at  each  j 
location  is  determined  by  the  relation 


r(y.u)2  zj*  + eD/2  (8) 

J  £  <Szuj«>2  +  E 

where  the  i  summation  is  carried  out  only  for  those  points  within  an  interval 

0.2D  <  z-  <  20.  Here  0  is  the  base  diameter,  6  a  central  difference  opera- 
J  '**  ^ 

tor,  and  e  a  positive  parameter  that  ensures  a  standard  grid  if  all  6Z u j ^  are 

zero  or  if  e  is  very  large.  Additional  averaging  is  used  in  the  x  direction 
(longitudinal  direction).  Each  of  these  grids  consisted  of  200  points  in  the 
longitudinal  direction  with  80  points  located  in  the  base  region  and  50  points 
in  the  normal  direction.  Details  of  the  grid  patching  used  consistent  with 
the  flow  field  segmentation  procedure  and  the  strategy  of  clustering  of  grid 
points  can  be  found  in  References  3  and  4. 


V.  RESULTS 


All  the  computations  were  made  at  Mm  =  2.0,  a  =  0  and  for  a  jet  exit 

Mach  number  of  2.5.  Solutions  were  obtained  for  the  jet-off  case  and  jet-on 
cases  for  jet-to-f ree-stream  pressure  ratios,  P . /P^  =  1,  3  and  9.  Figure  6 

J 

shows  a  schematic  illustration  of  the  base  flow  for  the  jet-on  condition  and 
its  associated  nomenclature. 

Figure  7  shows  the  velocity  vectors  in  the  base  region  for  the  jet-off 

case.  The  reci rculatory  flow  in  the  near  wake  is  clearly  evident.  The 

expected  velocity  defect  can  also  be  seen  further  downstream  in  this  figure. 

The  effect  of  the  centered  propulsive  jet  on  the  base  region  flow  field  is 

shown  in  Figures  8-10  for  P./P  =  1,  3,  and  9,  respectively.  Figure  8a  shows 

J 

the  velocity  vectors  in  the  near  wake  for  a  pressure  ratio  of  1.  The  flow 

field  in  the  base  region  has  changed  considerably  and  the  large  reci rcul atory 

bubble  seen  for  the  jet-off  condition  is  much  reduced.  This  small  separation 

bubble  is  a  region  of  counter  clockwise  recirculating  flow  and  can  be  seen 

near  the  base  corner  as  shown  in  Figure  8b.  Similar  features  can  be  seen  for 

the  pressure  ratio,  P./P  =  3  as  shown  in  Figure  9.  In  both  Figures  8  and  9, 

J  °° 

one  can  observe  the  oblique  compression  shock  wave  at  the  end  of  the  afterbody 
and  the  barrel  shock  inside  the  plume  indicated  by  the  turning  of  the  velocity 
vectors.  Figure  10a  shows  the  velocity  vectors  for  the  high  pressure  ratio 
case  of  P./P^  =  9.  The  shape  of  the  plume  is  clearly  shown.  An  expanded 

view  of  the  flow  field  near  the  end  of  the  afterbody  and  the  base  corner  is 
shown  in  Figure  10 b.  This  figure  shows  an  extensive  region  of  flow  separation 
upstream  of  the  base  corner.  The  small  separation  bubble  seen  downstream  of 
the  base  corner  for  the  lower  pressure  ratio  case  is  virtually  eliminated. 
The  separation  bubble  upstream  of  the  base  corner  is  confined  to  the  boundary 
layer  on  the  afterbody.  Additionally,  the  compression  shock  wave  seen  at  the 
end  of  the  afterbody  with  the  lower  jet  pressure  has  now  moved  further 
upstream  of  the  base  corner  with  the  high  jet  exit  pressure. 

The  next  four  Figures  11-14  show  the  comparisons  of  the  computed  density 
contours  and  the  experimental  Schlieren  pictures9  for  the  jet-off  case  and 
jet-on  cases  with  P./P  =  1,  3,  and  9,  respectively.  These  Schilieren  photo- 

j 

graphs  are  for  the  same  geometry  and  flow  conditions  obtained  from  the 
experimental  study  of  Agrell,  et  al.14  Figure  11  shows  the  comparison  for  the 
jet-off  condition.  For  this  case,  the  expansions  at  the  base  corner,  the  free 
shear  layer  and  the  recompression  shock  downstream  of  the  base  are  all  clearly 
observed  in  the  computed  results  and  agree  very  well  with  the  experimental 
observations.  Comparisons  for  the  jet-on  conditions  are  shown  in  Figures  12- 
14  and  it  is  clear  that  flow  features  in  the  base  region  have  changed  due  to 
the  presence  of  the  jet.  Figure  12  shows  the  results  for  the  jet  pressure 
ratio  of  1.  The  flow  features  to  be  seen  are  the  oblique  shock  at  the  end  of 
the  afterbody,  the  trailing  shock  system  inside  the  plume  and  the  slip  line 
that  emanates  from  the  nozzle  lip  and  defines  the  jet  boundary.  The  trailing 
shock  inside  the  plume  closes  about  1/2  caliber  downstream  of  the  exit  plane 
and  results  in  a  Mach  reflection.  As  the  jet-to-f ree-stream  pressure  ratio  is 
increased  to  3.0,  the  trailing  shocks  inside  the  plume  cross  each  other  about 
2  calibers  downstream  of  the  base  (see  Figure  13).  Other  features  are  similar 
to  the  pressure  ratio  of  one  case.  Details  of  the  flow  features  in  the  base 
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region  are  changed  as  the  jet  exit  pressure  ratio  is  increased  to  9.  The 
trailing  shock  system  seen  with  the  lower  pressure  ratios  is  not  observed  in 
either  the  computations  or  experiment  (see  Figure  14).  For  this  higher  jet 
pressure  ratio,  one  can  observe  a  lambda  shock  near  the  base  corner  which 
induces  a  separation  region  on  the  afterbody  just  upstream  of  the  base 
corner.  The  agreement  between  the  computed  and  the  experimentally  observed 
flow  features  is  very  good  for  both  jet-off  and  jet-on  conditions. 

Quantitatively,  one  is  interested  in  how  the  complex  flow  field  in  the 
base  region  affects  the  base  pressure.  Figure  15  shows  the  base  pressure  as  a 
function  of  boattail  angles  for  the  jet-off  condition.  Although  experimental 
measurements  were  made  for  various  boattail  angles,  computations  were 
restricted  to  only  the  8°  boattail  case.  The  computed  base  pressure  agrees 
well  with  the  experimental  data.  The  centered  jet  affects  the  base  region 
flow  field  considerably  and,  thus,  has  a  strong  affect  on  the  base  pressure. 
The  effect  of  the  jet  on  the  base  pressure  for  various  jet  to  free  stream 
pressure  ratios  is  shown  in  Figure  16.  The  trend  of  an  increase  in  base  pres¬ 
sure  with  jet  pressure  seen  experimentally  is  clearly  predicted  by  the  numer¬ 
ical  solutions.  Good  agreement  is  found  at  the  high  pressure  ratio.  Some 
disagreement,  however,  is  observed  at  the  lower  pressure  ratios. 


VI.  CONCLUDING  REMARKS 

A  computational  study  was  made  for  supersonic  flow  over  a  missile 
afterbody  in  the  presence  of  a  centered  propulsive  jet.  The  thin-layer  form 
of  the  compressible  Navier-Stokes  equations  was  solved  using  a  time-dependent, 
implicit  numerical  algorithm.  Solutions  were  obtained  for  both  jet-off  and 
jet-on  conditions  for  a  free  stream  Mach  number  of  2.0  and  the  jet  exit  Mach 
number  of  2.5.  Three  jet-to-f ree-stream  ratios  (1,  3,  and  9)  were 
considered.  The  grids  in  the  base  region  were  adapted  to  the  free  shear  layer 
as  the  solutions  developed.  Qualitative  features  of  the  base  region  flow 
field  such  as  the  compression  shock,  plume  shape,  and  trailing  shock  system 
seen  experimentally  were  easily  observed  in  the  computed  results. 
Quantitative  comparisons  indicate  good  agreement  for  the  jet-off  case.  Also, 
the  predicted  effect  of  the  jet  on  the  base  pressure  has  the  correct  trend 
observed  experimentally.  Some  disagreement  at  lower  pressure  ratios  exists. 
The  accuracy  of  these  predictions  should  improve  as  turbulence  modeling  used 
for  these  flows  improves. 
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Figure  7.  Velocity  vectors  in  the  base  region,  M  =  2.0,  jet-off. 


1.0 


Figure  12.  Computed  density  contours;  experimental  Schlieren 

photograph,  =  2.0,  Mj  =  2.5,  P./P^  =  1.0,  jet-on. 


1.3 


Figure  13.  Computational  contours;  experimental  Schlieren 

photograph,  =  2.0,  Mj  =  2.5,  P./Pot  =  3.0,  jet-c 


pressure  with  boattail  angle,  M  =  2.0,  jet-off 
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LIST  OF  SYMBOLS 


speed  of  sound 


cross  sectional  area  at  the  base 
base  drag  coefficient,  2  D^p^u^A 


specific  heat  at  constant  pressure 


pressure  coefficient,  2(p  -  Pj/p»u* 


body  diameter  (50mm) 
base  drag 


total  energy  per  unit  volume/p^a* 


E,  F,  q  flux  vector  of  transformed  Navier-Stokes 
equations 


n-invariant  source  vector 


Jacobian  of  transformati on 


Mach  number 
pressure/p  a* 


Prandtl  number,  u  c  /k 

OO  n  CD 


body  radius 

Reynolds  number,  P_,a  D/u 


viscous  flux  vector 


t  physical  time 

u,v,w  Cartesian  velocity  components/a 


LIST  OF  SYMBOLS 


U,V,W  contravari ant  velocity  components/a^ 

x,y,z  physical  Cartesian  coordinates 

a  angle  of  attack 

y  ratio  of  specific  heats 

<  coefficient  of  thermal  conductivity/^ 

u  coefficient  of  vi scosity/u^ 

5,n,c  transformed  coordinates  in  axial, 

circumferential  and  radial  directions 

p  density/p,,,, 

x  transformed  time 

<t>  circumferential  angle 

Superscri pt 

*  critical  value 

Subsen'  pt 
b  base 

j  jet  conditions 

J  longitudinal  direction 

L  normal  direction 

o  total  conditions 

st  stagnation  conditions 


ao 


free  stream  conditions 
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